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Abstract: We theoretically investigate the laser wavelength dependence 
of asymmetric dissociation of Hj. It is found that the electron localization 
in molecular dissociation is significantly manipulated by varying the 
wavelength of the driving field. Through creating a strong nuclear vibration 
in the laser-molecular interaction, our simulations demonstrate that the 
few-cycle mid-infrared pulse can effectively localize the electron at one 
of the dissociating nuclei with weak ionization. Moreover, we show that 
the observed phase-shift of the dissociation asymmetry is attributed to the 
different population transfers by the remaining fields after the internuclear 
distances reach the one-photon coupUng point. 

©2011 Optical Society of America 

OCIS codes: (020.2649) Strong field laser physics; (020.4180) Multiphoton processes; 
(190.4180) Multiphoton processes; (270.6620) Strong-field processes. 



References and links 

1. F. Krausz and M. Ivanov, "Attosecond physics," Rev. Mod. Phys. 81, 163 (2009). 

2. J. H. Posthumus, "The dynamics of small molecules in intense laser fields," Rep. Prog. Phys. 67, 623-665 (2004). 

3. A. D. Bandrauk, S. Chelkowski, and H. S. Nguyen, "Attosecond Localization of Electrons in Molecules," Int. J. 
Quant. Chem. 100, 834-844 (2004). 

4. V. Roudnev, B. D. Esry, and 1. Ben-Itzhak, "Controlling HD^ and dissociation with the canier-envelope 
phase difference of an intense ultrashort laser pulse," Phys. Rev. Lett. 93, 163601 (2004). 

5. M. F. Kling, Ch. Siedschlag, A. J. Verhoef, J. 1. Khan, M. Schultze, Th. Uphues, Y. Ni, M. Uiberacker, M. 
Drescher, F. Krausz, and M. J. J. Vrakking, "Control of electron localization in molecular dissociation," Science 
312, 246-248 (2006). 

6. X. M. Tong and C. D. Lin, "Dynamics of light-field control of molecular dissociation at the few-cycle limit," 
Phys. Rev. Lett. 98, 123002 (2007). 

7. S. Grafe and M. Yu. Ivanov, "Effective fields in laser-driven electron recollision and charge localization," Phys. 
Rev. Lett. 99, 163603 (2007). 

8. F. He, A. Becker, and U. Thumm, "Strong-field modulated diffraction effects in the correlated electron-nuclear 
motion in dissociating Hj," Phys. Rev. Lett. 101, 213002 (2008). 

9. M. Kremer, B. Fischer, B. Feuerstein, V. L. B. de Jesus, V. Sharma, C. Hofrichter, A. Rudenko. U. Thumm, C. 
D. Schroter, R. Moshammer, and J. Ullrich, "Electron localization in molecular fragmentation of H2 by carrier- 
envelope phase stabilized laser pulses," Phys. Rev. Lett. 103, 213003 (2009). 

10. D. Ray, F He, S. De, W. Cao, H. Mashiko, R Ranitovic, K. P Singh, 1. Znakovskaya, U. Thumm, G. G. Paulus, 
M. F. Kling, I. V. Litvinyuk, and C. L. Cocke, "Ion-energy dependence of asymmetric dissociation of D2 by a 
two-color laser field," Phys. Rev. Lett. 103, 223201 (2009). 

11. G. Sansone, F. Kelkensberg, J. F. Perez-Torres, F. morales, M. F. Kling, W. Siu, O. Ghafur, P. Johnsson, M. 
Swoboda, E. Benedetti, F. Ferrari, F. LePine, J. L. Sanz-Vicario, S. Zherebtsov, I. Znakovskaya, A. L'Huillier, M. 



Yu. Ivanov, M. Nisoli, F. Martin, and M. J. J. Vrakking, "Electron localization following attosecond molecular 
photoionization," Nature (London) 465, 763 (2010) 

12. K. P. Singh, F. He, P. Ranitovic, W. Cao, S. De, D. Ray, S. Chen, U. Thumm, A. Becker, M. M. Murnane, H. C. 
Kapteyn, I. V. Litvinyuk, and C. L. Cocke, "Control of electron localization in deuterium molecular ions using 
an attosecond pulse train and a many-cycle infrared pulse," Phys. Rev. Lett. 104, 023001 (2010). 

13. B. Fischer, M. Kremer, T. Pfeifer, B. Feuerstein, V. Sharma, U. Thumm, C. D. Schroter, R. Moshammer, and J. 
Ullrich, "Steering the electron in Hj by nuclear wave packet dynamics," Phys. Rev. Lett. 105, 223001 (2010). 

14. C. Calvert, R. B. King, W. A. Bryan, W. R. Newell, J. R McCann, J. B. Greenwood, and L D. Williams, "Multi- 
pulse scheme for enhancing electron localization through vibrational wavepacket manipulation," J. Phys. B 43, 
011001 (2010). 

15. F. Kelkensberg, G. Sansone, M. Y. Ivanov, and M. Vrakking, "A semi-classical model of attosecond electron 
localization in dissociative ionization of hydrogen," Phys. Chem. Chem. Phys. 13, 8647-8652 (2011). 

16. I. Znakovskaya, P. von den Hoff, S. Zherebtsov, A. Wirth, O. Herrweilh, M. J. J. Vrakking, R. de Vivie-Riedle, 
and M. F. Kling, "Attosecond control of electron dynamics in carbon monoxide," Phys. Rev. Lett. 103, 103002 
(2009). 

17. K. Liu, W. Hong, and P. Lu, "Phase dependence of electron localization in HeH^+ dissociation with an intense 
few-cycle laser pulse," Opt. Express 19, 20279 (201 1). 

18. C. Liu, Z. Zeng, P. Wei, P. Liu, R. Li, and Z. Xu, "Driving-laser wavelength dependence of high-order harmonic 
generation in Hj molecules," Phys. Rev. A 81, 033426 (2010). 

19. M. Magrakvelidze, F. He, T. Niederhausen, L V. Litvinyuk, and U. Thumm, "Quantum-beat imaging of the 
nuclear dynamics in D^: Dependence of bond softening and bond hardening on laser intersity, wavelength, and 
pulse duration," Phys. Rev. A 79, 033410 (2009). 

20. X. Gu, G. Marcus, Y. Deng, T. Metzger, C. Teisset, N. Ishii, T. Fuji, A. Baltuska, R. Butkus, V. Pervak, H. 
Ishizuki, T. Taira, T. Kobayashi, R. Kienberger, and F. Krausz, "Generation of carrier-envelope-phase-stable 2- 
cycle 740-^J pulses at 2.1-;tlm carrier wavelength," Opt. Express 17, 62 (2009). 

21. J. Tate, T. Auguste, H. G. Muller, P. Salieres, P. Agostini, and L. F. DiMauro, "Scaling of wave-packet dynamics 
in an intense midinfrared field," Phys. Rev. Lett. 98, 013901 (2007). 

22. Q. Zhang, P. Lu, P. Lan, W. Hong, and Z. Yang, "Multi-cycle laser-driven broadband supercontinuum with a 
modulated polarization gating," Opt. Express 16, 9795 (2008). 

23. C. Vozzi, F. Calegari, F. Frassetto, L. Poletto, G. Sansone, P. Villoresi, M. Nisoli, S. De Silvestri, and S. Stagira, 
"Coherent continuum generation above 100 eV driven by an ir parametric source in a two-color scheme," Phys. 
Rev. A 79, 033842 (2009). 

24. W. Hong, Q. Zhang, Z. Yang, and P. Lu, "Electron dynamic control for the quantum path in the midinfrared 
regime using a weak near-infrared pulse," Phys. Rev. A 80, 053407 (2009). 

25. K. C. Kulander, F. H. Mies, and K. J. Schafer, "Model for studies of laser-induced nonlinear processes in 
molecules," Phys. Rev. A 53, 2562 (1996). 

26. G. L. Ver Steeg, K. Bartschat, and I. Bray, "Time-dependent model calculations for a molecular hydrogen ion in 
a strong ultra-short laser pulse," J. Phys. B 36, 3325-3336 (2003). 

27. J. R. Hiskes, "Dissociation of molecular ions by electric and magnetic fields," Phys. Rev. 122, 1207-1217 (1961). 

28. B. Feuerstein and Uwe Thumm, "On the computation of momentum distributions within wavepacket propagation 
calculations," J. Phys. B 36, 707-716 (2003). 

29. V. Roudnev and B. D. Esry, "HD+ in a short strong laser pulse: Practical consideration of the observability of 
carrier-envelope phase effects," Phys. Rev. A 76, 023403 (2007). 



1. Introduction 

The interaction of intense laser pulses with atoms and molecules has been a very important 
topic in current research of attosecond physics. Recent advances in femtosecond (fs) laser tech- 
nology have made it possible for scientists to observe, study and control the ultrafast molecular 
dynamics in the strong field 1 1 1. Several mechanisms responsible for the laser induced dynam- 
ics of the molecule have been discovered in theoretical and experimental studies on H2 and its 
isotopes, such as bond softening, bond hardening, above threshold dissociation and so on Q. 

The electron localization in the dissociation of molecules has recently become an interest- 
ing subject for molecular research |j3]|4|. It is crucial for steering many physical and chemical 
processes to control the electron localization in molecules on the attosecond time-scale ||T). A 
landmark experimental demonstration in this direction was achieved by Kling et al. |j5), where 
a few-cycle carrier-envelope phase (CEP) stabilized pulse was used to control the localization 
of the electron in the dissociation of D^. As stimulated by the pioneer works |[3j-j5J, many 



methods have been proposed to control as well as explain the electron localization in small 
molecules |[6[|T5) and heteronuclear molecules | [T6][T7| . However, partly due to the limitation 
imposed by ultrafast amplifier technology, the wavelengths of the controlling pulse in those 
previous studies are mainly in the near-infrared region. Nonetheless, it has been shown that the 
dynamics of the nuclear wavepackets, which can be used to steer the localization of the electron 
in molecules fT3l, are significantly changed by mid-infrared pulses 1 18 19 1. It is therefore ap- 



pealing to know how the mid-infrared pulses influence the asymmetry of electron localization 
in molecular dissociation. Moreover, with the development of the infrared optical parametric 
amplifier system |20|, the mid-infrared pulses have been widely utilized to study as well as con- 
trol the electron dynamics in the ultrafast interaction |[2T[|24). The understanding of wavelength 
effect on the electron localization in molecular dissociation will provide helpful experimental 
guidance on steering physical and chemical processes in the future. 

In the present work, the dependence of the asymmetric dissociation of on the wave- 
length of the driving field is studied by numerical solutions of a non-Born-Oppenheimer time- 
dependent Schrodinger equation. The results show that both the amplitude and the relative 
phase of the CEP-sensitive asymmetry are changed significantly with the wavelength. We then 
qualitatively explain the wavelength-dependent dissociation asymmetry by analyzing the real- 
time propagation of the electron-nuclear wavepackets during the interaction. Further results 
demonstrate that the few-cycle mid-infrared pulse can create a strong nuclear vibration and 
effectively localize the electron at one of the dissociating nuclei with weak ionization. More- 
over, the observed phase-shift of the dissociation asymmetry is found to be attributed to the 
different population transfers by the remaining fields after the internuclear distances reach the 
one-photon coupling point. Although the recollision induced dissociation is demonstrated to be 
an important dissociation channel |5|^ and also expected to be wavelength dependent, it is not 
considered in the present study. 

2. Theoretical model 

For numerical simulation we have solved the non-Born-Oppenheimer time-dependent 
Schrodinger equation (TDSE) for a reduced-dimensionality model of Hj. The model consists 
of one-dimensional motion of the nuclei and one-dimensional motion of the electron ||25][26|. 
In this model, the electronic and nuclear motions are restricted along the polarization direction 
of the Unearly polarized laser pulse. The time-dependent Schrodinger equation is written as 
(Hartree atomic units are used throughout unless otherwise indicated): 

= [Ho + V{t)]^{R,z;t), (1) 

where Hq is the field-free Hamiltonian, 

I I 1 1 1 

and y(f) is the electric potential including the laser-molecular interaction. Here, R is the in- 
ternuclear distance, z is the electron position measured from the center-of-mass of the protons, 
jj. = mp/2 and fi^ — 2mp/{2mp + 1) are the reduced masses with nip the mass of the pro- 
ton, and a — I and /3 = 0.03 are the soft-core parameters. The interaction with the laser is 
given by (in the dipole approximation and length gauge) | [27) V{t) = [1 + {2mp + i)^^]zE{t) 
with E{t) = £'oexp[— 21n(2)(f/T)^] sin{(Ot + (j)). Eq is the peak electric field amplitude, T is the 
pulse duration, (O is the central frequency, and (j) is the CEP of the pulse. The electric field is 
chosen to be linearly polarized along the z axis. Then one can solve Eq.(l) by using the Crank- 
Nicolson method. Compared to the 3D TDSE methods used in Refs. ||4]|8][T7l, the electric 




Fig. 1. (a) The probability density of tlie Hj molecular ground state and the configuration 
space regions associated with different final state channels. The dash lines indicate the 
positions of the nuclei. The color scale is linear, (b) The sketch of the potential energy 
surface for the IjcTg and 2pcT„ states of H^. 



motion is restricted along the z direction for the present approach, and consequently the elec- 
tronic wavepackets spread in the radial direction (perpendicular to the z-axis) is not taken into 
account. But it has been shown before that this type of reduced-dimensional model reproduces 
experimental result at least qualitatively |25 26) . 



In the present simulation, the wave function is discrete on a two-dimensional grid which 
ranges in /^-direction from to 32 a.u., for z from -300 to 300 a.u. with 800 and 3000 points, 
respectively. We have defined different regions of the grid to distinguish the dissociation, ion- 
ization and bound probabilities. As illustrated in Fig.[TJa), the bound region is defined as 

Pb : < < 10 and - 10 < z < 10, (3) 

and the two channels of dissociation are defined as 



P+:R>\Q and ^{z-R/2)^<5, (4) 
P : R>m and J{z + R/2)^ < 5, (5) 



where P+ and P indicate the probabilities of the electron localization on the upper and lower 
protons (p), respectively, and thus these are the directions that the H atom is emitted. By inte- 
grating the density of the final wavepackets in the regions described above, the final localization 
and bound probabilities can be obtained. The region desciibing the ionization is not contained 
in Fig. [T]but we estimate the ionization probability by storing the absorbed wavepackets reach- 
ing the boundary in the z direction. In our calculation, the ground state of (electronic IsGg 
state and vibrational v = state), as shown in Fig. [Tf a), is chosen to be the initial state of the 
system and obtained by propagating the field-free Schrodinger equation in imaginary time. The 
evolution of the wave function in the external field has been continued with a time step of 
dt = 0.05 until the probability of the bound region is converged. During the propagation, we 
have employed absorbing boundaries using cos '/^-masking functions but stored the absorbed 
contributions as dissociation or ionization, respectively. 

For better understanding of our simulation, here we briefly introduce the origin of the electron 
localization asymmetry in molecular dissociation. Figure [TJb) shows the sketch of the potential 
energy surface for the IsGg and 2pa„ states of H^- These two states are the gerade and ungerade 
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Fig. 2. The dependence of the phase-sensitive asymmetry parameter A on the wavelengths 
of 800, 1000, 1200 and 1600 nm (in increasing order from the left to the right column) for 
three different intensities of 8 x lO'^ (top row), 6 x lO''* (middle row) and 4 x lO'^ W/cm^ 
(bottom row). 

States, respectively, and the superposition of the dissociative wavepackets on the two states will 
lead to the electron localization on one of the nuclei. Depending on the CEP fT|, intensity ||8) 
and even wavelength of the pulse, population is transferred between the two states, changing 
the electron localization between the nuclei. The detail of the physical mechanism responsible 
for the asymmetry of electron localization will be discussed in Section 3. 

3. Results and discussion 

We aim to study the dependence of the asymmetric dissociation of on the wavelength with 
an intense two-cycle (FWHM) laser pulse. Let us first define the absolute asymmetry of the 
dissociation as A = f + — f _ p4) . Here, the absolute asymmetry parameter A is utilized because 
the total ionization and dissociation probabilities are comparable in our simulation and we at- 
tempt to find out the maximum probabilities of both the dissociation and electron localization. 
Therefore, in the present work this parameter A gives a more transparent representation of the 
asymmetric dissociation, whereas a normalized parameter could provide a misleading indica- 
tion; for example, a large normalized asymmetry could be obtained in a scenario where the total 
dissociation probability is negligible. 

Figure |2] reveals the dependence of the CEP-sensitive asymmetry parameter A (diamonds) 
on the wavelength for three intensities of 8 x 10^'*, 6 x 10^^ and 4 x 10^^ W/cm^, respectively. 
On the one hand, the amplitude of the asymmetry is changed significantly with the wavelength. 
For higher intensity (top row), the maximum dissociation asymmetry A„,av decreases with the 
wavelength, while for lower intensity (bottom row), we observe an opposite variation of A,„ax- 
Here, A„,„x is defined as the maximum of |A| in each panel. For the intermediate intensity 
(middle row), A^ax first increases and then decreases with the wavelength. On the other hand, 
the data of the asymmetry parameter A in Fig. [2] are fitted to sine functions (thick curves). We 
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Fig. 3. (Multimedia online) Snapshots of electron-nuclear probability density distribu- 
tion taken at f = 30 fs for four cases, in which the laser parameters for (a) (Media 1), (b) 
(Media 2), (c) (Media 3) and (d) (Media 4) are the same as those in Figs.|2fa),|2];d),|2i) 
and[2jl), respectively, except that the CEP is fixed at 0.5;r. The color scale is logarithmic. 



find that there is a phase-shift of the sine function between the panels. This means that A^ax 
appears at different values of the CEP when the wavelength or intensity is changed. 

In order to gain an intuitive understanding of the electron localization in the dissociation of 
HJ, we simulate the real-time propagation of the electron-nuclear wavepackets in the external 
fields of four cases, as shown in Figs.|3ja) to[3|d) (Media 1-^), respectively. For the four cases, 
the CEP of the pulse is fixed at 0.5;r, and the other laser parameters are the same as those in 
Figs.|2]^a),|2|d),|2|i) and|2|l), respectively. By following the evolution of the probability density 
in the external fields (Media 1^), a strong vibration of the probability density can be observed 
when the pulse is applied. Meanwhile, part of wavepackets are ionized in the ±z direction at 
the extrema of the electric field. Apparently, the ionized probability density for higher intensity 
is much more than that for lower intensity. At the tail of the pulse, the ionization almost stops. 
Pait of the vibrational wavepackets dissociate slowly along the upper or lower proton, and the 
rest stay in the bound states. 

According to the real-time evolution of wavepackets in the interaction, the dependence of 
Amax on the wavelength for the relative higher and lower intensities in Fig.|2]can be qualitatively 
explained as follows, hi the case of 8 x lO'^ W/cm^, a great deal of the wavepackets are ionized. 
For 800 nm, most of the wavepackets dissociate with the protons due to the nuclear vibration 
after the ionization [Fig. |3ja). Media 1]. But longer wavelength will cause stronger vibration 
and thus much more wavepackets will be further ionized by the driving field of high intensity 
pS) . Consequently, for 1600 nm, almost all the wavepackets have been ionized before the 
pulse is turned off and only few wavepackets are dissociated [Fig.|3jb), Media 2]. Note that the 
absolute asymmetry parameter A is closely associated with the total dissociation probability 
| [l4) . Thus, under the intensity of 8 x lO'^ W/cm^, the reduce of the dissociation probability 
leads to the decrease of A,„ax with the wavelength. However, the situation is quite different in 
the case of relative lower pulse intensity, where few wavepackets are ionized due to the relative 
low intensity. For 800 nm, because of the weak nuclear vibration, there are few dissociative 
wavepackets and most of the wavepackets stay in the bound states [Fig. [3]^c), Media 3]. But 




Fig. 4. Time-dependent ionization and dissociation probabilities under four wave- 
lengths for peak intensities of 4 x lO'^ (left column), 6 x lO''* (middle column) and 
8 X lO'"* W/cm^ (right column). The color arrows indicate the dependence of the final 
probabilities on the pulse wavelength. 



the nuclei vibrate more and more intensely with increasing wavelength. For 1600 nm, only a 
few wavepackets remain in the bound states. Most of them are dissociated and localized on the 
dissociating protons [Fig.[3]^d), Media 4]. As a result, A,„aj( increases with the wavelength for 
the intensity of 4 x lO'"^ W/cm^. 

To gain further insight in the relation between the ionization and dissociation that influences 
the maximum asymmetry of electron localization, the time-dependent ionization and dissocia- 
tion probabilities during the interactions are calculated. The results under four wavelengths for 
different intensities are depicted in Fig.|4] Firstly, in the left column of Fig.|4] the total ioniza- 
tion probabilities for all wavelengths are less than 15%, and the total dissociation probability 
increases rapidly with the wavelength. However, when the intensity increases, the ionization 
increases more rapidly for longer wavelength, as shown in the middle column of Fig. |4] More 
than 40% of the wavepackets are ionized for 1200 nm and 1600 nm, which directly reduce the 
total dissociation probabilities. Lastly, the ionization are further enhanced for the intensity of 
8 X lO'^ W/cm^, with the result that the total dissociation probability decreases with wave- 
length, as shown in the right column of Fig. |4] Note again that the asymmetry A is closely 
related to the dissociation probability. Thus the variation of the total dissociation probability 
under different intensities discussed above leads to the corresponding dependence of Amux on 
the wavelength in each row of Fig. [2] 

Now we discuss the wavelength effect on the electron localization in more detail by compar- 
ing the results in Fig.|2]and Fig.|4] On the one hand, as shown in Fig.|4|e), the final dissociation 
probabilities for 800 nm (solid curve) and 1200 nm (dotted curve) are approximately equal. 
However, the corresponding A„,av for 1200 nm in Fig. |2|g) is larger than that for 800 nm in 
Fig. |2je). On the other hand, although the dissociation probability for the wavelength of 800 
nm under the intensity of 4 x lO'"* W/cm^ [solid curve in Fig.Uld)] is larger than that for 1600 
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Fig. 5. The final dissociation (dots), ionization (stars) and survival (rings) probabilities as a 
function of the intensity for the wavelengths of (a) 800, (b) 1000, (c) 1200 and (d) 1600 nm. 
The CEP is fixed at 0.5;i:. 

nm under 8 x lO"^ W/cm^ [dash-dot curve in Fig. ^f)], the A„k„ in the two cases are almost 
the same, as shown in Figs.|2];i) and|2|d). In addition, some similar situations can be also found 
for other laser parameters. These features show that the longer-wavelength pulses are more 
effective to create stronger electron localization asymmetry. 

For comprehensive understanding of the laser-molecular interaction under different pulse 
wavelengths, we show the dissociation, ionization and bound probabilities as a function of the 
pulse intensity varying from 5 x 10^^ to 1 x 10^^ W/cm^ for the wavelengths of 800, 1000, 
1200 and 1600 nm in Figs.|5ja) to|5jd), respectively. The CEP of the pulse is fixed at Q.5n and 
the results do not vary significantly with the CEP in our calculation. Obviously, there is a peak 
value of the dissociation probability for each wavelength in Fig. [5] and its position moves from 
high to low intensities when the wavelength is increased. The offset of the dissociation peak 
can be understood as that the longer- wavelength pulse can create stronger nuclear vibration and 
thus the dissociation probability increases more rapidly. But meanwhile, because of the intense 
nuclear vibration for the long-wavelength pulse, the ionization also increases rapidly, resulting 
in the early decrease of the dissociation probability. 

By focusing on the given region of the intensity between the vertical dash lines in Fig. |5] 
one can find that the variations of the dissociation probability with the intensity under differ- 
ent wavelength are generally coincident with those of A„,av in each column of Fig. |2j respec- 
tively. For the laser parameters of Figs. |2|a), |2|f) and|2]^l), strong asymmetries are obtained, 
approximately corresponding to the dissociation peaks in Fig. |5] However, when the dissoci- 
ation reaches its peak value, the ionization probability is quite high for 800 nm and becomes 
lower for longer wavelength. For the case of 1600 nm and 4 x lO'^ W/cm^ in Fig.jsjd), the total 
dissociation and ionization probabilities are 0.60 and 0.14, respectively. By varying the CEP of 
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Fig. 6. CEP-dependent and CEP-averaged KER spectra for different combinations of the 
pulse wavelength and intensity. The laser parameters for (a), (c) and (e) are the same as 
those in Figs.|2|a), (f) and (1), respectively. 



the pulse, the A,„ax shown in Fig.|2|l) are as high as 0.273, where the electron localization prob- 
abilities P+ and P are 0.436 and 0.163, respectively. In this scenario, more than 72% of the 
dissociative (and 43% of all) wavepackets are localized on one of the dissociating nuclei. This 
indicates that the mid-infrared pulse can effectively localize the electron at a chosen nucleus 
with weak ionization. 

Next, we turn to the physical mechanism responsible for the phase-shift of the asymmetry 
A in Fig. [2] First of all, to obtain more detailed insight into the CEP control of electron local- 
ization, we calculate the fragment kinetic energy release (KER) distribution using the 'virtual 
detector' method p8] . Figures [6| a), [6|c) and|6|e) show the CEP-dependent KER distributions 
for different combinations of the pulse wavelength and intensity, respectively. The solid curves 
in the Figs.|6|b),|6|d) and|6|f) indicate the corresponding CEP-averaged KER spectra. As shown 
in Fig.|6j the KER distributions peak at around 2 eV, 1 .4 eV and 0.6 eV for the 800-nm, 1000-nm 
and 1600-nm pulses, respectively, indicating that the above-threshold dissociation contributes 
significantly to the KER spectra. The contributions from bond-softening are relatively weak 
due to the high pulse intensity (> 4 x 10^'* W/cm^) |29) . 

According to the KER spectra in Fig. |6j we propose a possible simplified explanation for 
the phase-shift of the asymmetry. As illustrated in Fig. [TJb), the wavepackets are bound in 
the initial state at first and begin to vibrate after the pulse is applied. Then some vibrational 
wavepackets on the IsGg state absorb three photons and transfer to the dissociative IpGu state. 
When the wavepackets on both states pass through the one-photon crossing point, the electric 
field couples the two states and then population is transferred between them. The coherent 
superposition of the dissociative wavepackets on the gerade and ungerade states finally leads 
to a localization of the electron at either of the nuclei |8 -10|. Because (a) the motion of the 
wavepackets depends on the pulse intensity and wavelength, (b) the one-photon crossing point 
depends on the wavelength, and (c) the coupling field is sensitive to the CEP, the population 
transfer will be different if the intensity, wavelength or CEP is changed. As a result, the electron 
localization depends on the CEP, intensity and wavelength of the pulse. 

To verify the interpretation for the origin of the phase-shift discussing above, we show a 
comparison of the time-dependent asymmetry versus the time-dependent internuclear distance 
for different laser parameters in Fig.|7] Here are the calculation details. For the TDSE approach 
beyond Born-Oppenheimer approximation in the present study, the individual time-dependent 
population on the ground or excited states couldn't be calculated in principle. Instead, we can 
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Fig. 7. The electric fields, internuclear distances and asymmetries as a function of time (in 
different units for the left and right columns). For the left column, the pulse parameters are 
with the same wavelength of 1600 nm and intensity of 4 x lO'^ W/cm^ but different values 
of the CEP. For the right column, the pulse parameters are with the same CEP of 0.5;r but 
different wavelengths and intensities. The horizontal and vertical dash lines indicate the 
one-photon coupling points (about 6.1 a.u. for 1600 nm and 5.2 a.u. for 800 nm) and the 
time when the average internuclear distances reach the points, respectively. 



extract the time-dependent bound wavepackets by 

1'bound(/?,z;f) - L MR,z)\1'{R,z;t))v„{R,z), (6) 

,1=0 

with v„ the vibrational bound states and = 21 the number of v„ in our model. Then the time- 
dependent probabilities of the dissociative wavepackets localized on the upper or lower nucleus 
can be estimated by 

rlO r±lQ 

P'^{t) ^ dR dz|^(/;,z;f)-^bound(-R,z;f)l' 







dR dz|^(7?,z;f)-^bound(-R,z;f)h (7) 

10 J±f-5 

where R,„ax corresponds to the boundary of the grid. The time-dependent expectation values of 
the internuclear distance R{t) for the total wavepackets is given by 1 26 1 

J,'-dRj-'^2^dzRmR,z-,t)\^ 
!o^-dRj^-^dz\^{R,z;t)\^ 

with ±Zmax the boundaries of the grid. 

We firstly focus on the situation for different values of CEP and otherwise identical laser 
parameters, as shown in the left column of Fig.|7] The variations of the time-dependent internu- 
clear distances are almost coincident before they reach the coupling point, and the asymmetries 



[P+{t) — Pi (f)] oscillate with the electric fields in a similar way. However, after the internuclear 
distances reach the coupling point, the asymmetries vary in completely different ways and fi- 
nally lead to the different localizations of the electron, as shown in Fig.[7]^c). Then we turn to 
the situation for the same CEP but different wavelengths and intensities, as shown the right 
column of Fig.|7] The optical cycle (o.c.) has been used as the unit of time and the intensities of 
the pulses have been normalized so that the electric fields share the same profile. In Fig. |7|e), 
the internuclear distances reach their coupling points approximately at f = 2 o.c. (but actually 
not the same time). Before t = 2 o.c, both the asymmetries in Fig.[7|f) oscillate with the pulse 
almost in a same way. However, because of the different molecular motions as well as pho- 
ton energies, the asymmetries change quite differently after the coupling point, resulting in the 
electron localization on different nuclei. Overall, the two situations discussing above indicate 
that the one-photon coupling point plays an important role in determining the asymmetry, and 
the asymmetry depends not only on the CEP but also the intensity and wavelength of the pulse 
even though the electric field share the same profile. 

4. Conclusion 

In conclusion, we have studied the wavelength dependence of the asymmetric dissociation of 
HJ by solving numerically the non-Born-Oppenheimer time-dependent Schrodinger equation. 
By varying the wavelength, significant changes of the electron localization asymmetry have 
been observed. The strong asymmetries that are obtained at different intensities are associated 
with the offset of the dissociation peak, and the phase-shift of the asymmetry is attributed to the 
different population transfers by the remaining fields after the internuclear distances reach the 
one-photon coupling point. Moreover, a high degree of control over the electron localization on 
a chosen nucleus with relative low ionization probability of the molecule has been achieved with 
the 1600-nm pulse, suggesting that the mid-infrared few-cycle pulse can be used to effectively 
control the electron localization in molecular dissociation with weak ionization. 
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